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Recent  geological  and  geodetic  studies  have  suggested  that  the  region  suirounding  Vandenberg  AFB  is  undergoing 
active  crustal  deformation,  with  important  implications  for  both  the  geodetic  stability  and  the  seismogenic  pmential 
of  the  Western  Test  Range  (WTR)  {Feigl  ei  a!..  1990).  Part  of  the  evidettcc  for  significant  deformation  was  obtained 
from  GPS  measurements  over  a  broad  area  of  central  and  southern  California,  which  we  carried  out  in  cooperation 
with  other  university  and  government  scientists  from  1987  through  1991.  Although  useful  in  defining  the  regional 
tectonic  setting,  these  measurements  are  of  insufficient  spatial  and  temporal  density  to  answer  many  imptxtant 
questions  about  the  seismogenic  potential  of  Vandenberg.  The  sites  observed  in  the  March  1992  survey  are  given  in 
Tabte  I,  and  a  map  of  the  enlarged  network  i.s  shown  in  Figure  1.  Four  new  sites  (CASM,  FAARF,  RDRK,  and 
FIGP)  were  est^Ushed  on  or  bridging  three  of  tlie  four  major  anticlines  that  cut  the  Santa  Maria  Basin.  Our  i»imary 
scie*itific  goal  is  to  localize  the  measured  deformation  on  one  or  more  of  these  structures.  One  of  the  new  sites 
(SOAP)  iHOvides  a  mme  stable  anchor  for  that  part  of  our  netwok  south  of  the  Santa  Ynez  River  Fault  Three  other 
sites  occupied  for  the  first  time  (ARG3.  VINA.  VANP),  as  well  as  two  previously  occupied  (ALVA,  VNDN),  arc 
existing  DMA  sites  near  the  south  Base  PGG  A  site  (VNDP).  This  small  sub-network  provides  a  means  to  moniUH'  the 
local  stability  of  the  three  sites  (VAND,  VNDN.  and  VNDP)  that  have  been  or  will  be  used  for  long-term  monitoring 
of  defixmation.  The  two-day  occupation  of  the  South  Base  sub-network  also  ixovided  us  the  tqiportunity  to  study 
more  carefully  the  effects  of  atmospheric  water  '^por  on  GPS  measurements. 


98  5)  ,31  o5 


.»jr-rTr  .-.rr.-.-  —  ii mi r^rmt-g^xTTirr-nr.T- 


OTRiMtr 

(U) 


09  mt  MOi 

(U) 


Ml  uaflTAiiwQMiiriucr 


Crustal  Deformation  Measurements  in  the  Vicinity  of 
Vandenberg  Air  Force  Base 

Grant  AFOSR-90.0339 
(MIT  OSP  No.  75390) 

Annual  Technical  Report 
for  the  period 

1  September  1991  -  14  October  1992 

Submitted  to 

Air  Force  Office  of  Scientific  Research 

Stanley  Dickinson 
Program  Manager 

AFOSR/NL  Bolling  AFB,  DC  20332-6448 

Vernita  J.  Slater 
Contracting  Officer 
AFOSR/PKZA 

Bolling  AFB,  DC  20332-6448 


QTTAUnr  INSPECTED  4 


Robert  W.  King 
Principal  Investigator 

Department  of  Earth,  Atmospheric,  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 


f 


2 


Introduction 

Recent  geological  and  geodetic  studies  have  suggested  that  the  region  surrounding 
Vandenberg  AFB  is  undergoing  active  crustal  deformation,  with  important  implications  for 
both  the  geodetic  stability  and  the  seismogenic  potential  of  the  Western  Test  Range  (WTR) 
[Feigl  et  al,  1990].  Part  of  the  evidence  for  significant  deformation  was  obtained  from 
GPS  measurements  over  a  broad  area  of  central  and  southern  California,  which  we  carried 
out  in  cooperation  with  other  university  and  government  scientists  from  1987  through 
1991.  Although  useful  in  defining  the  regional  tectonic  setting,  these  measurements  arc  of 
insufficient  spatial  and  temporal  density  to  answer  many  important  questions  about  the 
seismogenic  potential  of  Vandenberg. 

In  1989  we  received  funding  from  AFOSR  under  the  Commander’s  Reserve  Fund 
(Grant  AFOSR-89-04(X)),  with  matching  funds  from  MIT,  to  purchase  five  GPS  receivers 
and  to  begin  a  series  of  measurements  designed  to  determine  the  laagnitude  and  spatial 
distribution  of  deformation  in  the  Santa  Maria  Fold  and  Thrust  Belt  (SMFTB),  a  region 
encompassing  the  major  faults  and  folds  within  50  km  of  Vandenberg.  Two  additional 
receivers  were  purchased  in  May  1992  and  have  since  been  installed  in  continuously 
operating  GPS  stations  at  Vandenberg  and  the  China  Lake  Naval  Weapons  Center.  They 
are  operated  as  part  of  the  Permanent  GPS  Geodetic  Array  (PGGA)  in  California, 
providing  the  ability  to  monitor  not  only  interseismic  deformation,  but  any  transient  motion 
which  might  occur  prior  to,  during,  or  after  an  earthquake. 

The  pre-grant  studies  of  Vandenberg  tectonics  and  our  progress  during  the  first  year  of 
this  grant  have  been  completely  described  in  the  doctoral  thesis  of  Kurt  L.  Feigl,  completed 
at  M.  I.  T.  in  September  1991  and  submitted  with  last  year’s  annual  report. 

GPS  Measurements 

Since  Febmary  1990,  we  have  carried  out  five  GPS  surveys  involving  the  stations  of 
the  Vandenberg  network.  Three  of  these  surveys  (February  and  September  1990,  and 
March  1992)  included  most  or  all  of  the  current  network;  the  other  two  (March  1990  and 
February  1991)  included  two  or  three  Vandenberg  stations  as  part  of  a  remeasurement  of 
the  regional  central  and  southern  California  networks.  The  details  of  the  first  four  surveys 
are  given  in  Chapter  4  of  Feigl  [1991],  especially  Figures  4.1-4.3,  Tables  4.1-4.6,  and  the 
text  on  pp.  124-125  and  137.  In  the  fifth  survey  (March  1992)  we  remeasured  the  relative 
positions  of  eight  of  the  Vandenberg-network  stations  surveyed  in  1990,  and  also 
established  nine  new  stations  to  densify  the  network. 

The  sites  observed  in  the  March  1992  survey  are  given  in  Table  1 ,  and  a  map  of  the 
enlarged  network  is  shown  in  Figure  1.  Four  new  sites  (CASM,  FARF,  RDRK,  and 
FIGP)  were  established  on  or  bridging  three  of  the  four  major  anticlines  that  cut  the  Santa 
Maria  Basin.  Our  primary  scientific  goal  is  to  localize  the  measured  deformation  on  one  or 
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more  of  these  structures.  One  of  the  new  sites  (SOAP)  provides  a  more  stable  anchor  for 
that  part  of  our  network  south  of  the  Santa  Ynez  River  Fault.  Three  other  sites  occupied 
for  the  first  time  (ARG3,  VINA,  VANP),  as  well  as  two  previously  occupied  (ALVA, 
VNDN),  are  existing  DMA  sites  near  the  South  Base  PGGA  site  (VNDP).  This  small 
subnetwork  provides  a  means  to  monitor  the  local  stability  of  the  three  sites  (VAND, 
VNDN,  and  VNDP)  that  have  been  or  will  be  used  for  long-term  monitoring  of 
deformation. 

The  two-day  occupation  of  the  South  Base  subnetwork  also  provided  us  the 
opportunity  to  study  more  carefully  the  effects  of  atmosperic  water  vapor  on  GPS 
measurements.  Vandenberg  is  a  promising  setting  for  this  study  since  it  is  subjected  to 
strong  atmospheric  gradients  from  the  onshore  winds.  Both  very  long  baseline 
interferometry  (VLBI)  and  GPS  observations  conducted  there  over  the  past  decade  show 
atmospheric  errors  that  are  greater  than  those  for  many  other  sites  in  Calfomia  sites.  Our 
study  is  further  enhanced  by  the  availablility  of  radiosonde  data  obtained  twice  daily  by  the 
base  weather  detachment.  Although  the  detachment  normally  alternates  launches  between 
North  and  South  Vandenberg,  during  the  period  of  our  observations  they  altered  their 
schedule  to  take  all  the  measurements  from  South  base. 


Table  1.  Observations  acquired  during  the  March  1992  survey 


Receiver 


Mar 

day 

day 

MITl 

MIT2 

MIT3 

MIT4 

CITl 

CIT2 

Mon 

2 

062 

ALAM 

- - 

— 

— 

VAND 

VNDN 

Tue 

3 

063 

AlAM 

— 

GRAS 

LIND 

VAND 

VNDN 

Wed 

4 

064 

LOSP 

— 

GRAS 

LIND 

VNDP 

VNDN 

Thu 

5 

065 

LOSP 

— 

GRAS 

LIND 

VNDPa 

VNDNa 

Fri 

6 

066 

LOSP 

— 

FIGP 

— 

- a 

- a 

Sat 

7 

067 

FARF 

DMAV 

FIGP 

MADC 

VNDPa 

- a 

Sun 

8 

068 

FARF 

VNDN 

VINA 

MADC 

VNDP 

SO.AP 

Mon 

9 

069 

RUSl 

3ARG 

VINA 

ALVA 

VNDP 

SOAP 

Tue 

10 

070 

RUSl 

3ARG 

VINA 

ALVA 

VNDP 

SOAP 

Wed 

11 

071 

— 

— 

VINA 

— 

VNDP 

— 

Thu 

12 

072 

LIND 

DMAV 

RDRK 

— 

VNDP 

VNDN 

Fri 

13 

073 

LIND 

DMAV 

P^DRK 

CASH 

VNDP 

VNDN 

Sat 

14 

074 

LIND 

ALAM 

RDRK 

CASH 

VNDP 

GRAS 

Sun 

tot 

15 

075 

Litro 

13 

AlLAM 

8 

12 

CASM 

12 

VNDP 

10 

GRAS 

10  64 

2a  1 

Days  069  and  070  constitute  the  "atmospheric  experiment", 
a  =  cow  attack 


Figure  1.  Tectonic  map  of  Vandenberg  and  vicinity  (Santa  Maria  Fold  and  Thrust  Belt)  showing  the 
location  of  faults  (thick  lines)  and  folds  (thin  lines)  and  of  the  GPS  stations  (dots  with  4 -letter  codes) 
included  in  our  study.  Solid  dots  indicate  stations  which  have  been  surveyed  both  in  1990  and  1992;  open 
dots  indicate  stations  surveyed  for  the  first  time  in  1992.  Only  three  of  the  seven  stations  in  the  dense 
South  Base  network  near  Point  Arguello  are  shown. 
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Data  Analysis 

GPS  observations  in  the  vicinity  of  Vandenberg  cannot  be  understood  in  isolation  from 
the  rest  of  southern  California.  Only  by  determining  (vector)  motions  of  the  stations  in  the 
Vandenberg  network  relative  to  stations  outside  the  network  can  we  assess  the  relative 
importance  of  dilatation  and  rotation  in  the  pattern  of  deformation.  Moreover,  it  is  essential 
to  include  in  the  analysis  observations  from  a  continental-  or  global-scale  networic  in  order 
to  reduce  the  errors  due  to  unmodeled  motions  of  the  GPS  satellites.  Thus,  we  have 
analyzed  all  of  the  Vandenberg  data  simultaneously  with  both  GPS  and  \1.BI  data  from 
central  and  southern  California  and  a  global  network  of  stations.  Much  of  our  effort  during 
the  past  year  has  been  spent  in  refining  this  analysis  and  describing  it  in  a  comprehensive 
paper  submitted  in  December  to  the  Journal  of  Geophysical  Research  [Feigl  et  ai,  1992]. 
Figure  2  shows  our  estimates  of  the  velocities  of  35  sites  in  central  and  southern  C^ifomia. 
The  dominant  motion  at  this  scale  is  the  simple  shear  due  to  the  San  Andreas  Fault  (SAF) 
system,  which  contributes  a  few  millimeters  per  year  even  to  sites  as  far  west  as 
Vandenberg. 

In  order  to  understand  better  more  local  sources  of  deformation,  we  show  in  Figure  3 
residual  velocities  after  removing  an  approximate  model  of  velocities  associated  with  the 
SAF.  Of  the  five  stations  in  the  Santa  Maria  Fold  and  Thrust  Belt  (SMFTB),  Mount  Lospe 
(LOSP)  and  Madre  (MADC)  exhibit  residual  motions  significantly  different  from  zero  at 
95%  confidence.  In  particular,  the  residual  velocity  of  Madre  with  respect  to  the 
Vandenberg  South  Base  site  (VNDN)  implies  1.8  ±  0.9  mm/yr  of  shortening  and  2.9  ±  0.9 
mmy/yr  of  right-lateral  slip  in  the  frame  defined  by  the  local  (N60“W)  strike  of  the  folds. 

The  rate  of  shortening  is  compatible  with  the  geological  rate  of  2-5  mm/yr  estimated 
from  a  balanced  cross  section  [Namson  and  Davis,  1990],  but  smaller  than  the  6  ±  1  mm/yr 
estimated  from  our  earlier  comparison  of  GPS  and  historical  survey  data  [Feigl  et  ai, 
1990].  The  latter  study,  however,  assumed  uniform  strain  in  the  SMFTB,  no  net  rotation, 
and  no  strain  accumulating  from  the  locked  SAF.  All  three  of  these  hypo'Jieses  are  suspect 
in  light  of  the  spatially  varying  strain  rates  and  the  suggestion  of  rotation  from  our  analysis 
of  the  larger  California  network,  and  of  our  model  calculation  of  3  mm/yr  of  relative  motion 
between  Vandenberg  and  Madre. 

The  amount  of  right-lateral  strike-slip  shear  in  the  SMFTB  is  larger  than  expected, 
given  the  lack  of  faults  active  in  the  Quaternary  [Jennings,  1975;  Sylvester  and  Darrow, 
1979;  Clark  et  ai,  1984].  One  explanation  would  be  unmodeled  strain  accumulation  on  the 
SAF.  Alternatively,  strain  accumulation  on  the  offshore  Hosgri  fault  [Hall,  1987;  1981] 
may  be  indicated.  The  addition  of  data  from  the  March,  1992,  survey  and  the  continuous 
observations  from  the  Vandenberg  PGGA  site  should  help  to  distinguish  among  these 
hypotheses. 


Figure  3.  Residual  vlocities  with  respect  to  Vandcnberg  (VNDN)  after  removal  of  a  model  of  shear  strain 
due  to  the  San  Andreas  fault  system,  from  FeigI  et  al.  (1992J.  Note  that  the  velocity  Kale  is  double  that 
of  the  previous  figure.  The  tectonic  features  include  the  San  Andreas  fault  (SAF),  the  Santa  Barbara 
Channel  (SBC),  and  the  southern  Coast  Ranges  (SCR). 
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An  unexpected  result  from  the  Vandenberg  PGGA  station  has  been  the  contribution  of 
its  observations  to  geodetic  determination  of  far-field  displacements  firwn  the  Landers  (Mw 
7.3)  and  Big  Bear  (Mw  6.2)  earthquakes  of  28  June  1992.  These  earthquakes  occured 
within  100  km  of  the  PGGA  sites  at  Pinyon  Flat  Observatory  (PINl),  operated  by  the 
University  of  C^ifoma  at  San  Diego,  and  the  NASA  Goldstone  Complex  (GOLD)  (Figure 
4).  These  two  sites,  plus  the  PGGA  sites  at  Scripps  (SI02)  and  the  Jet  Propulsion 
Laboratory  in  Pasadena  (JPLM),  were  displaced  by  more  than  a  centimeter  by  the 
earthquakes.  Our  analysis  of  these  data,  in  collaboration  with  our  colleagues  at  Scripps, 
has  been  reported  in  a  letter  to  Nature  published  this  month  and  attached  to  this  report 
[Bock  et  aL,  1993].  A  supplemental  plot  of  the  Vandenberg  displacements,  not  included  in 
the  Nature  letter,  is  shown  in  Figure  5.  There  is  a  hint  of  a  displacement  of  a  few 
millimeters  in  both  the  north  and  east  component  of  the  Vandenberg  position.  The  north 
displacement,  estimated  formally  to  be  4.5  ±  1.6  mm,  is  significantly  larger  than  the  value 
(0.9  mm)  predicted  by  our  model.  An  analysis  of  the  Vandenberg  observations  over  a 
longer  period  of  time  will  be  necessary  to  determine  if  the  uncertainty  in  our  estimated 
displacement  is  realistic.  Note  also  in  Figure  5  the  large  (>50  mm)  systematic  signature  in 
the  vertical  position  of  Vandenberg.  We  believe  that  this  is  due  to  deficiencies  in  our 
models  for  satellite-orbit  and  atmospheric  effects.  The  long  series  of  PGGA  observatons 
provide  an  important  data  set  for  improving  these  models. 


Figure  4.  Observed  (solid  arrows)  and 
modelled  (blank  arrows)  displacements  of 
PGGA  stations  duing  the  Landers/Big  Bear 
earthquakes,  shown  with  95%  confiducencc 
ellipses.  The  contours  are  for  the 
displacement  magnitude  based  on  an  elastic 
half-space  model  fit  to  measured  near-ficld 
displacements  and  the  inversion  of  strong- 
motion  seismic  data.  Also  shown  are  the 
estimates  obtained  from  an  independent 
analysis  of  the  PGGA  data  by  Blewitt  et  al. 
[1993]  at  JPL.  See  Bock  et  al.  [1993]  for 
details. 
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Figure  5.  Time  series  of  daily  positions  of  the 
Vandenberg  PGGA  antenna  with  respect  to  a  global 
reference  frame.  The  error  bars  reivesent  our  best 
estimate  of  the  one-standard  deviation  uncertainty  of 
each  estimate  and  have  been  scaled  to  account  for 
correlated  noise  sources.  This  series  may  be 
compared  directly  with  those  given  in  Figure  3  of 
Bock  et  a}.  {1993]  for  other  sites  in  the  PGGA. 
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LETTERS  TO  NATURE 


Few  of  the  trends  are  significant,  large  surface  warmrng  trends 
are  not  observed,  and  like  the  Western  Arctic  Ocean  analysis 
in  Table  1.  significant  surface  cooling  trends  are  found  during 
winter  and  autumn.  Significant  warming  trends  are  observed  at 
the  850  hPa  level  and  the  850-700 hPa  layer  during  winter,  in 
agreement  with  GCM  simulations,  but  the  surface  trend  is 
negative.  The  trends  presented  in  Table  2  are  more  representative 
temporally,  as  they  use  36-38  years  out  of  a  possible  41,  as 
compared  with  20-30  years  in  the  regional  analysis  (Table  1). 
On  the  basis  of  these  tests,  we  feel  that  any  possible  bias 
introduced  by  the  non-uniform  database  is  small. 

The  lack  of  widespread  significant  warming  trends  leads  us 


Receiveii  23  July,  jccepteo  2r  Oecemtier  1992 

1  HougMon.  J  t  IM  *  Otmaie  Ctange  Trie  PCC  Soenlific  Assessmeni  (Cambudge  umv  Press 
Canwridge.  I990i 

2  BuUyko.M  I  Telius  11.  611-619(19691 

3,  Seders.  W  D  2  epol  MeteOfOl  •,  392 -dOO  119691, 

4  mgrarri.  W  i  Wilson,  C  A  *  Milchell.  F  J  B  2  geiyytys  «es  94,  8609  8622  (19891 

5  Ke»y  P  M  Jones  P  D  Sear.  C  B .  Cherry.  B  S  G  4  TavePol.  8  K  Mon  Weam  110,  71  83 
(19821 

6  Hansen.  J  4  Lebeoeff.  S  2  geaptiys  Res  »2,  13345-13372  •.987| 

7  Hansen,  j  er  *  J  gecphys  Res  *3,  9341  -9364  (1988; 

8  Wa(sh  2  E  4  Crane  R  G  Geophya  Res  Lett  1».  29  32  (1992i 

9  Wa(sh.  J  E  4  Chaoman,  W  L  2  Climete  i,  237-250  (19901 
10  *n*e((.  J  K  Man  wealft  Rev  Hi,  1922  1930  (19861 

n  AnjeK  J  H  4  Horshover.  J  Mon  Weith.  Rev  111.  901-921  ';983l 
13  Raroly.  D  J  Geopnys  Res  Lett  16.  465  468  (1989J 

13  Kara.  JO  et  91  J  gecphys  Res  (submitleiJ) 

14  Iimerev  A  4  £  Egorov  s  A  UelerorOI  OiOre/  No  7  50  56  •.9911 


to  conclude  that  there  is  no  strong  evidence  to  support  mode! 
simulations  of  greenhouse  warming  over  the  Arctic  Ocean  for 
the  period  1950-1990.  Our  results,  combined  with  the  incon¬ 
sistent  performance  of  model  simulations  of  Arctic  climate'’, 
indicate  a  need  to  understand  better  the  physical  processes  that 
affect  the  polar  regions,  especially  atmosphere  ice-ocean  inter 
actions,  ocean  heat  transfer  and  cloud  radiative  effects,  and  to 
incorporate  thermodynamic  sea  ice  components  into  future 
models.  We  consider  the  retrieval  of  temperafure  profiles  from 
satellite  sounding  instruments  to  be  an  important  means  of 
further  resolving  the  spatial  and  temporal  gradients  in  Arctic 
air  temperatures.  □ 


15  A  P,  Ttmerev,  A  A  4  feorov.  S  A  AJ^aa  Sam  SSSfi  Dont  mo  1U3 

16  SwTc/c.  M  C .  Kahl,  j  0  4  Schneif  R  C  J  Oimaie.  8.  6 J  5  629  (1 99?) 

1?  Setreie.M  C  et  at  J  geophys  Res  97.  9411  9422  (1*192) 

18  J  6  Atoft  Wests  ffev  108.  1527-1535  11977) 

19  »«.JM  )  D .  Set^e^e.  M  C .  S .  Skooy  S  m  4  Scnrwfi,  C  BuW  4m  mef  Soc  73. 

1824  1830  <1992) 

20  Seftete  M  C .  KaN.  J.  0  4  Shotani  S  Sstro^  mkJ  tee  Oats  Ce-'^tef  spec  Rep  Sc  2 
'CtRES.  UmvefSity  of  Colorado.  1992) 

21  S  M  thesis.  Untv  o<  Wisconsm-Mtiwaiitee  {1992> 

22  Oaconis.  P  4  Efron.  B  Scte*>t  Am  248.  116-130(1983) 

AC*<sCWlEDGEACfvTS  Portions  of  this  wod*  we'O  soonsored  Oy  the  hOAA  s  Oirna'.e  and  O^obai 
0<r4«  Program  toe  Eiecirtc  Power  Research  )ns-*!ute  and  (he  D>vis«on  of  Poiar  Prt^rarns 
Dfc.ec*  (teas  cortducted  as  part  of  the  US-Rossia  Committee  on  Cooperation  n  Env>ronmenta« 
P'C'ec'on.  the  influence  of  Environmentaf  Change  or  Climate  (Wortrmg  Group  Vlli)  Hetpfui  comments 
e  p'OviOed  Oy  A  A  Tsonis 


Detection  of  crustal  deformation 
from  the  Landers  earthquake 
sequence  using  continuous 
geodetic  measurements 

Yehuda  Bocit,  Duncan  C.  Agnew,  Peng  Fang, 

Joachim  F.  Genrich,  Bradford  H.  Hager*, 

Thomas  A.  Herring*,  Kenneth  W.  Hudnutt, 

Robert  W.  King*,  Shawn  Larsent,  J.-Bernard  Minster, 
Keith  Stark,  Shimon  Wdowinski  &  Frank  K.  Wyatt 

institute  of  GeopHysics  and  Planetary  Physics,  Scr  ops  Institution  of 
Oceanography.  La  Joila.  California  92093.  USA 
‘Department  o'  Earth,  Atmospheric,  and  Planeta'.  Sciences. 
Massachusetts  I'-stitute  of  Technology.  Cambridge 
Massachusetts  02139.  USA 

tUS  Geological  Survey,  Pasadena,  California  91 106  USA 
T Scientific  SoftAare  Division,  Lawrence  Livermore  National  Laboratory 
Livermore.  California  94550.  USA 


The  measurement  of  crustai  motions  in  tectonically  active  regions 
is  being  performed  increasingly  by  the  satellite-based  Global 
Positioning  System  (CPS)’'^  which  offers  considerable  advantages 
over  conventional  geodetic  techniques^'*.  Continuously  operating 
GPS  arrays  with  ground-based  receivers  spaced  tens  of  kilometres 
apart  have  been  established  in  central  Japan’*  and  southern 
California  to  monitor  the  spatial  and  temporal  details  of  crustal 
deformation.  Here  we  report  the  first  measurements  for  a  major 
earthquake  by  a  continuously  operating  CPS  network,  the  Per¬ 
manent  GPS  Geodetic  Array  (PCGA)’"’  in  southern  t.'alifnrnia. 
The  Landers  (magnitude  Af„  of  7.3)  and  Big  Bear  (Af,  6.2) 
earthquakes  of  28  June  1992  were  monitored  by  daily  observations. 
Ten  weeks  of  measurements,  centred  on  the  earthquake  events, 
indicate  significant  coseismic  motion  at  all  PCGA  sites,  significant 
post-seismic  motion  a:  one  site  for  two  weeks  after  the  earthquakes, 
and  no  significant  preseismic  motion.  These  measurements  demon¬ 
strate  the  potential  of  CPS  monitoring  for  precise  detection  of 
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precursory  and  aftershock  seismic  deformation  in  the  near  and 
far  field. 

The  PGGA,  established  in  southern  California  (Fig,  1)  in  the 
spring  of  1990,  is  a  network  of  five  continuously  opera'ing  GPS 
receivers  providing  an  unim  rrrupted  record  of  crustal  motion 
in  near  real-time.  At  each  site  there  is  a  precise  P-code  GPS 
receiver  with  its  antenna  mounted  on  a  geodetic  monument. 
Twenty-four  hours  of  data  are  automatically  collected  from  each 
site  once  a  day;  the  operational  analysis  provides  the  site  position 
averaged  over  the  day  (0-24  h  Universal  Time  Coordinated, 
UT(T.  although  finer  resolution  is  possible.  The  precision  of 
the  daily  relauve  position  between  any  two  sites  in  the  network, 
based  on  long-term  scatter  of  nearly  2  years  of  measurements, 
i,  ~5mm  in  the  horizontal  and  10-20  mm  in  the  vertical.  To 
obtain  this  precision,  and  to  achieve  near-real-time  solutions, 
we  compute  orbits  for  all  GPS  satellites,  using  data  collected 
by  a  globally  distributed  network  of  about  25  permanent  tracking 
stations'"  (Fig.  I),  and  corrections  to  tabu'ated  predictions"  of 
the  orientation  of  the  Earth's  rotation  axis  (polar  motion).  The 
worldwide  tracking  network  defines  a  global  reference  frame  in 
which  coordinates  for  the  PGGA  stations,  in  unstable  southern 
C  alifornia,  can  be  computed  with  respect  to  rigid  plate  interiors 
(for  example,  the  North  American  plate). 

Here  we  report  measurements  of  seismically  induced  displace¬ 
ments  of  the  PGGA  stations  due  to  the  Landers  {M^  7,3,  1  i ;  58 
UTC,  34.22°  N,  116.43°  W)  and  Big  Bear  (Af,.  6.2,  15  :07  UTC, 
?4  21°  N,  1 16.83°  W)  earthquakes  of  28  June  1992.  We  examined 
the  series  of  PGGA  station  positions  in  the  iO-week  period 
centred  on  the  day  of  the  earthquakes  using  a  Kalman-filter 
formulation'-  to  analyse  the  daily  24-hour  PGGA  solutions.  For 
the  day  of  the  earthquakes,  we  computed  the  station  positions 
separately  from  the  12  hours  of  data  before  the  Landers 
earthquake  and  the  9  hours  of  data  after  the  Big  Bear  event 
The  global  tracking  network  was  fairly  extensive  as  the  3-month 
Imernational  T-PS  Service  campaign"  began  on  21  June  (Fig. 
1 1.  Displacements  of  the  PGGA  sites,  listed  in  Tabic  I,  were 
determined  by  examining  the  variation  in  the  positions  before 
and  after  the  earthquakes  with  respect  to  a  global  reference 
frame  defined  by  the  coordinates  of  the  tracking  network 
stations.  The  horizontal  displacements  of  the  PGGA  sites  arc 
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FIG  1  The  dislribution  of  glooai  GPS  permanent  tracking 
stations  and  PG6A  sites  used  (or  our  analysis  of  observations 
for  the  period  May-August  1992  The  closest  tracking  site 
outside  California  is  DfiAO  m  Penticton.  British  Columoia  All 
stations  use  Rogue  SNR -8  GPS  receivers  except  for  two  PGGA 
sites.  S102  in  la  Jolla  and  VNOP  at  Vandonoerg  An  Force 
Base,  which  use  Ashtech  iX  il3  GPS  receivers  Trie  global 
ti  acking  network  is  usually  described  by  the  acronyms  ClGNET 
and  FIMM  standing  for  Cooperative  international  GPS  Network 
and  Fiducial  Latxjratones  (or  an  International  Natural  science 
Network'®,  and  more  recently  by  the  International  GPS  Service 
(IGS)" 


IGS  tracking  network 

■y 


plotted  in  Fig.  2.  The  maximum  horizontal  displacement  of 
48  mm  was  detected  at  PINS,  located  at  the  Pinon  Flat  Observa¬ 
tory  (PFO)  ~80  km  from  the  seismic  rupture  zone.  In  Fig  3  we 
take  differences  for  the  positions  obtained  with  respect  to  the 
global  reference  frame  and  plot  the  daily  record  of  relative 
horizontal  positions  between  PI  N 1  and  three  other  stations  (Fig 
1).  These  were  a  global  tracking  site  (DRAG)  near  Penticton, 
British  Columbia,  Canada,  more  than  1,700  km  to  the  north 
and  the  global  tracking  station  closest  to  the  PGGA;  the  PGGA 
site  (JPLM)  in  Pasadena;  and  the  PGGA  site  (GOLD)  at  the 
NASA  Deep  Space  Network  Goldstone  Complex. 

An  examination  of  the  station  displacements  indicates  that 
most  of  the  deformation  is  due  to  the  coseismic  phase  of  the 
crustal  deformation  cycle'''  The  coseismic  displacements  appear 
clearly  as  step  functions  in  the  time  series  of  daily  station 
positions  ( Fig.  3 ).  No  significant  pre  -seismic  signature  is  discern¬ 
ible  from  the  five  weeks  of  daily  data  taken  before  the 
earthquakes.  There  appears,  however,  to  be  a  significant  posi- 
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seismic  signature  of  0  9*0  3  mm  day  ‘  in  our  estimates  of  the 
relative  positions  of  GOLD  and  PIN  I  for  16  days  after  the 
earthquakes  (Fig.  3)  An  examination  of  the  displacements  at 
the  individual  stations  indicates  that  most  of  the  post-seismic 
motion  occurs  at  GOLD  at  a  rate  of  0.7  ±0.3  mm  day  (We 
were  concerned  that  displacements  observed  at  GOLD  mighi 
stem  from  the  construction  of  its  antenna  support,  a  12-m  high 
microwave  tower.  This  station,  unlike  the  other  PGOA  stations 
which  have  very  stable  geodetic  monuments,  was  neither  instal¬ 
led  nor  operated  primarily  for  monitoring  tectonic  motions,  and 
might  easily  be  recovering  from  the  mainshock  accelerations 
Comparing  the  changing  position  of  this  monument  w  ith  respect 
to  another  continuously  operating  GPS  receiver  10  km  away 
indicates,  however,  that  the  displacements  were  tectonic  in 
origin  1  The  remainder  of  the  motion,  '•3  mm  in  total,  occurs 
a;  PINl,  but  the  rate  of  displacement  is  not  well  determined 
from  the  data.  Greater  temporal  detail  and  higher  shon-term 
tesolution  is  provided  by  laser  strainmeters  at  PFO;  these  show 
r.o  pre-seismic  signal,  aiiliougii  post-scismic  s:  in  the  first 

1-2  weeks  following  the  earthquakes  are  consist. nt  with  our 
obvervations  at  PI  Nl.  Funher  discrimination  of  the  post  sei>mic 
signal  will  require  the  implementation  of  more  refined  analysis 
techniques  We  intend  to  study  the  much  longer  time  series  of 
data  before  and  after  the  earthquakes  to  look  for  possible 
pre-seismic  signals  and  changes  in  the  interseismic  rate  of  defor- 
-jtion.  and  to  understand  better  the  error  spectrum  of  con- 
t  nuous  GPS  data  We  arc  currently  investigating  an  increase  in 


-32  Observed  (solid  arrows)  and  modelled  (blank  arrows)  displace "ents 
i'  Ike  PGGA  stations  mcluamg  95%  confiderice  ellipses  The  observed 
■:  soiarements  are  calculated  with  respect  to  the  reference  frame  oe'ned 
c.  the  positions  and  velocities  o*  the  global  tracking  stations,  and  me  ode 
fe  total  displacement  estimated  over  the  5-week  period  attc'  tne 
ea  thquakes  Except  for  GOLD  and  a  very  small  effect  at  PINl,  the  disc  ace- 
~ents  seem  to  be  due  entirely  to  coseismic  motion  Tne  total  displacement 
ccserved  at  GOLD  is  17  mm  which  ncludes  -6  mm  of  coseismic  rnotion 
Z'd  Hmm  of  apparent  post-seismic  motion  over  the  16  days  following 
earthquakes.  For  comparison  we  stvow  the  dispiacements  and  95% 
ccnfidence  ellipses  (unshaded)  computed  by  Blewitt  et  al  '^.  The  contours 
c'  displacement  magnitude,  and  the  calculated  displacements  are  fc'  an 
e  astic  halfspace  (all  units  are  millimetres)  The  surface  trace  of  the  U'ders 
■■p'ute  (heavy  line)  is  composed  of  six  segments  with  end  points  at  lat'tjde 
8"d  longitude  (degiees)  (34  150)0.243  57081.  (34  2250  243  5€?5i 
34  3625  243  53751,  (34  4 5ai  243  5000),  (34  4958  243  4-50' 

34  5917, 243  3917)  and  (34  6542  243  3208)  The  magnitudes  Of  s''C  are 
'.5  3  5, 1 0.  4  8.  5  2  and  0  8m  right-lateral  respectively  with  all  segments 
'•cm  0  10  15  km  depth.  The  dashed  line  shows  the  surface  trace  assumed 
'C'  the  fault  segment  of  the  Big  Bear  earthquake  with  0  4  m  of  left  la'e  al 
s  c  on  8  segment  from  3  ;o  15km  depth  and  end  cents 
34  1200. 243  1037)  and  (34  2884  233  3297) 
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LETTERS  TO  NATURE 


^ABLE  1  Tneofeiical  ami  oDserved  displacements 


Site 


PSN1 

Gao 

JPtM 

SI02 

VNOP 


Distance 

(km) 

89 

117 

155 

185 

380 


North  (mm) 

East  (mm) 

AmpiituOe  (mm) 

0 

M 

0 

M 

0 

42  8 

456!  12 

16  2 

140  t  21 

45  8 

47  7  z  24 

-175 

-146  1 1  4 

-8  5 

-8  0  1  2  4 

19  5 

l66z  28 

4  7 

34tl3 

-13  4 

-14  7  ±24 

142 

151 t  2  7 

151 

130i I 7 

60 

101 ±28 

162 

165x3  3 

09 

45716 

-30 

-41 ±25 

31 

6  1  ±  3  0 

Anmuth 

Veftcal  (rnm. 

WO 

(si 

0 

M  0 

M 

0 

0  96 

20  7 

17  1 

36 

138 

96  !  64 

117 

205  9 

208  7 

2  8 

7  3 

60i  56 

0  94 

289  3 

2830 

63 

0 

0 

0  98 

21  7 

376 

161 

89 

10  3  !  8  2 

0  51 

286  7 

317  7 

-310 

26 

39x54 

Total  displacements  in  the  hotizontal  components  of  the  PGGA  stations  computed  by  a  variable  slip  dislocation  model  (M)  and  estimated  horn  10  weeks 
of  GPSobseivations  (0)  centred  on  the  day  of  the  earthquakes,  and  a  comparison  of  the  modelled  and  observed  amjMiludes  and  a/imuins  of  the  dtspiacements 
The  displacements  are  computed  with  respect  to  the  global  reference  frame  defined  by  the  global  trackirg  stations  (Fig  i).  The  errors  listed  tw  the  observed 
displacements  are  1  <7  values  The  distance  from  the  PGGA  sites  to  the  centroid  of  the  model  is  given  For  completeness  we  include  the  vertical  tisplacements 
computed  from  the  dislocation  model  and  estimated  from  the  (5PS  observations  in  this  case,  the  ladte  gives  the  v^ues  for  the  motion  of  the  PGGA  sues 
relative  to  JPIM  which  has  been  shifted  to  0  for  both  observation  and  model  values  This  requires  addir^  3  7  mm  to  the  model  value  and  25,9  mm  to  tre 
observed  value  The  vertical  component  errors  are  not  well  understood,  however,  so  the  agreement  with  the  dislocation  model  may  be  Irxtuitous 


the  data  noise  after  the  earthquakes  (Fig.  3i.  this  may  partly  be 
a  meteorological  effect  related  to  the  transition  into  summer,  a 
phenomenon  that  has  been  observed  in  many  years  of  positions 
determined  by  very-long-baseline  interferometry  in  California 
An  independent  analysis  of  the  PGGA  and  global  tracking 
data  by  a  group  at  the  Jet  Propulsion  Laboratory  (JPL)  yielded 

DRAG  to  PIN  1  (Length  1 ,738  989  km) 


GOLD  to  PIN1  (Length  L05  004  km) 


Day  number,  1992 


FIG.  3  Time  series  of  daily  horizontal  relative  positions  (north  and  east 
components)  over  a  10- week  period  centred  on  the  day  of  the  earthrjuakes 
(day  180,  28  June  1992)  for  PINl  relative  to  OSAO  (a,  b).  JPLM  Ic.  d)  and 
Gao  (e.  i)  The  dotted  horizontal  lines  are  determined  from  the  weighted 
means  of  the  data  for  the  5-week  period  before  ard  after  the  earthquakes, 
and  indicate  the  coseismic  signature  (these  differ  from  the  Kalman-filter 
estimates  given  in  Table  1  by  rx>t  more  than  2.5  mm).  The  error  bars  have 
been  scaled  by  a  factor  of  two  to  account  for  the  effects  of  correlated  ntjise 
sources  These  are  not  included  in  our  statistical  esf  mates  of  the  uncertaint¬ 
ies,  which  assume  white-noise  error  sources  In  e  we  fit  by  ieasl  squares 
a  straight  line,  with  a  sitjpe  of  0,9 ±0  3 mm  day  v  to  the  apparent  post- 
seismic  signature  for  the  16  days  following  the  earthquakes  The  east 
component  has  a  larger  scatter  than  the  north  component  because  the 
Integer -cycle  phase  ambiguities  were  not  resrjived  to  integer  values'*'*  in 
the  daily  P<3GA  solutions 


similar  total  displacements,  as  reponcd  by  Blewitt  ei  of".  The 
JPL  group  used  the  GIPSY  GPS  software'*”  to  process  the 
data  whereas  we  used  the  GAMIT/CI.OBK  GPS  software''". 
The  algorithms  used  in  these  two  programs  were  indeprendently 
developed  and  use  different  approaches.  Furthermore,  the  two 
groups  used  different  global  reference  frames  and  different  sub¬ 
sets  of  the  global  tracking  and  PGGA  data.  We  analysed  ten 
weeks  of  data  from  all  five  PGGA  sites,  which  included  three 
Rogue  SNR-8  GPS  receivers  and  two  Ashtech  LX-II3  GPS 
receivers  (at  VNDP  and  S102I,  whereas  Blewitt  et  at.  analysed 
8  w  eeks  of  Rogue  receiver  data  only.  A  comparison  of  the  results 
is  shown  in  Fig.  2.  The  overlaps  of  our  respective  95%  confidence 
ellipses  indicate  good  agreement  at  the  three  common  sites. 

In  Fig.  2  we  show  the  coseismic  displacements  computed  from 
a  variable  slip  dislocation  model  in  an  elastic  halfspace”-*’.  The 
model  for  the  Landers  earthquake  separates  the  6S-km  rupture 
into  six  venical  planar  segments  (Fig.  2),  oriented  to  coincide 
with  the  epicentre  and  the  curvilinear  pattern  observed  in  the 
surface  break  and  aftershock  distribution''.  The  right-lateral  slip 
on  the  southernmost  three  segments  was  calculated  by  averaging 
over  each  segment  the  observed  surface  offset  and  the  slip 
distribution  obtained  by  inverting  strong-motion  seismic  data''. 
Slip  on  the  northern  three  segments  was  obtained  by  fitting  the 
displacement  vector  at  GOLD.  This  required  significantly  less 
slip  than  observed  along  the  surface  break.  The  surface  faulting 
observed  geologically  continued  10  km  to  the  northwest  of  the 
main  aftershock  zone,  suggesting  that  rupture  along  this  segment 
of  the  fault  is  confined  to  shallow  depths'*''.  The  Big  Bear 
event  was  modelled  as  left-lateral  rupture  along  a  vertically 
dipping  fault,  oriented  to  coincide  with  the  focal  mechanism 
and  aftershock  distribution;  the  slip  was  chosen  to  agree  with 
the  seismic  moment  (H.  Kanamori,  personal  communication). 

There  is  good  agreement  between  the  theoretical  and  observed 
displacements  at  the  PGG  \  sites  (Tzble  !,  Fig.  2).  Our  fit  to 
the  displacement  vectors,  although  nonunique,  suggests  less  slip 
in  the  northern  half  of  the  Landers  rupture.  In  fact,  the  geodetic 
moment  calculated  using  these  far-field  PGGA  data  is  0.8  x 
IO'’'’Nm,  which  is  generally  less  than  moments  obtained  by 
independent  means;  near-field  seismic  (0. 8-0.9  x  lO'®  N  m),geo 
logical  (0.9  X  I0’“  N  m),  near-field  geodetic  (l.Ox  jO'®  N  m)  and 
teleseismic  (1.1  x  lO'®  N  m)''  '’.  The  comparatively  low  PGG/. 
moment  may  be  a  manifestation  of  the  influence  of  the  mantle, 
which  has  higher  values  of  elastic  moduli  than  the  crust",  so 
that  far-field  amplitudes  of  theoretically  computed  displace 
mems  are  reduced  relative  to  homogeneous  halfspace  models 
In  this  case,  the  far-field  is  distances  greater  than  the  30  km 
thickness  of  the  crust. 

These  data  show  that  continuous  GPS  arrays  can  provide 
reliable,  precise  and  rapid  determination  of  crustal  motion,  in 
panicular  seismically  induced  deformation.  Although  dense 
spatial  coverage  with  such  stations  is  not  economically  feasible 
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On  28  June  1992,  the  largest  earthquake  in  ralitornia  in  40  years 
(surface-wave  magnitude  Af.  =  7.5)  occurred  near  the  small  town 
of  Landers,  in  southeastern  California,  and  was  rollowred  three 
hours  later  by  the  nearby  Af.  6.5  Big  Bear  earthquake'.  For¬ 
tuitously,  the  Landers  earthquake  sequence  coincided  with  the  first 
week  of  the  official  three-month  test  period  of  the  International 
Global  Positioning  System  and  Geodynamics  Service^  (IGS). 
giving  us  an  unprecendented  opportunity  to  detect  absolute  pre-. 
CO-  and  post-seismic  displacements  at  a  distance  of  50-200  km 
from  the  main  rupture  with  millimetre-level  precision.  Mutual  and 
independent  confirmation  of  some  of  our  geodetic  results  are 
demonstrated  by  Bock  er  al.  in  this  issue'.  For  the  Landers 
earthquake,  the  observed  displacements  indicate  that  the  depth  of 
the  bottom  of  the  rupture  is  shallower  towards  the  northern  end, 
displacements  were  dominantly  symmetric,  and  the  rupture  exten¬ 
ded  further  south  on  the  Johnson  Valley  fault  than  has  been 
mapped  on  the  basis  of  surface  ground  offsets.  The  combined 
geodetic  moment  for  the  Landers  and  Big  Bear  earthquakes  (LI  x 
10“*  N  m"^')  agrees  well  with  teicseismic  estimates. 

The  Landers  and  Big  Bear  earthquakes  and  their  aftershocks 
occurred  along  faults  that  form  a  triangle  bounded  to  the  south¬ 
west  by  the  San  Andreas  fault  ( Fig.  1 ).  F.stensive  surface  rupture 
resulting  from  these  events  has  been  reported  along  the  Johnson 


Valley  and  Camp  Rock/tmerson  faults  Ground  breakage 
occurs  along  a  70-km  stretch  of  these  faults,  reaching  a  man.-num 
surface  offset  of  6.7  m  (ref  U  The  I  anders  earthquake  oveurred 
toward  the  southern  end  of  the  hypothesized  'Mojave  shear 
zone’,  which  trends  N.IS^W  across  the  Mojave  Desen,  inio 
Owens  Valley  and  the  northern  Basin  and  Range  province*  ’. 
This  zone  reportedly  carries  7-8  mm  yr  '  of  the  telatne  rriotion 
between  the  Pacific  and  North  American  plates,  and  ma,  be  a 
manifestation  of  a  suberusu!  fault*  '.  Aftershocks  following  the 
Landers  earthquake  tine  up  along  this  apparent  shear  zone  from 
as  far  south  as  the  San  Andreas  fault,  to  further  north  than  our 
station  GOLD  shown  in  Fig.  1.  The  pattern  of  aftershocks  is 
sparse  on  the  Camp  Rock  fault  segment  which  underlies  the 
northernmost  13  km  of  the  visible  surface  rupture*.  Using  a  new 
geodetic  tool  for  earthquake  studies,  we  have  estimated  per¬ 
manent  surface  displacements  in  southern  California  due  to  the 
cumulative  effect  of  events  on  28  June,  and  show  that  geodetic 
methods  provide  valuable  information  on  aspects  of  the  rupture 
mechanism  not  available  with  other  techniques.  We  also  offer 
an  explanation  for  the  unexpected  lack  of  aftershocks  on  the 
Camp  Rock  fault. 

Since  21  June  1992,  a  worldwide  network  of  stations  has  been 
routinely  receiving  precise  microwave  ranging  data  from  the 
United  States  Department  of  Defence's  IS-satellitc  Global 
Positioning  System  (GPS)’,  and  transmitting  the  data  to  IGS 
data  centres  to  be  made  available  to  analysis  centres  and  geo¬ 
dynamics  researchers.  Regional  GPS  networks  benefit  from  the 
precise  orbit  determination  and  reference  frame  stability  sup¬ 
plied  by  an  extensive  tracking  network*"*.  A  regional  array  of 
receivers  operated  jointly  by  the  Jet  Propulsion  Laboratory  (JPL) 
and  Scripps  Institution  of  Oceanography  has  been  operational 
in  southern  California  since  1990  (ref.  11).  A  simultaneous 
analysis  of  GPS  data  from  the  California  array  combined  with 
the  global  network  data  has  allowed  us  to  estimate  the  absolute 
displacements,  in  the  international  terrestrial  reference  frame'^ 
ilTRF),  of  three  stations  located  within  50-200 km  from  the 
Landers  earthquake  rupture,  with  2  mm  precision  in  the 
horizontal  plane. 

To  reduce  systematic  errors  that  can  Oe  introduced  by  mixing 
different  types  of  GPS  receivers  and  antennas,  we  have  analysed 
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"Kj  1  Map  showing  absolute  motions  ot  Goldstore  (GOLD).  Pasadena  (JflM). 
Pinyon  Flat  (PIN!)  and  Deadman  (DEAD)  Solid  arrows  are  the  observed 
displacements  with  95%  confidence  regions  Tbe  vectors  and  confidence 
region  for  DEAD  is  shown  at  01  limes  the  scale  of  the  other  stations  Tne 
model  displacements,  assuming  an  elastic  half  space,  are  shown  as  dashed 
arrows  The  surface  trace  of  the  model  of  the  LSnders  earthquake  is  she  -n 
by  the  solid  heavy  line  Dashed  heavy  line  (BB)  is  the  trace  of  the  fault  used 
to  model  the  Big  Bear  earthquake  Shaded  solid  and  dashed  lines  are  active 
'auits  in  the  region  JV  is  the  Johnson  Valley  fault:  CRE  is  the  Camp 
Rock/Emerson  fault 
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